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FOREWORD .. V . 

This report was prepared by the MRD Division of the General American 
Transportation Corporation under Contract OCD-OS- 63^176 for the Office of 
Civil Defense. Mr. Frank C. Allen of OCD's Directorate of Research was 
project monitor. This report 'covers the period of June 1963 to August 1964. 

The contract was a study of the. habitability of typical identified shelter 
configurations with consideration given to various space- allotments, ventilating 
rates, environmental criteria,land structural particulars.' All of,these parum- 
eters were studied over a wide range of variation with the primary emphasis on 
aboveground shelter spaces. 


Analytical studies of the environment in shelters are continuing under 
a subcontract with Stanford Research Institute, and a comprehensive report 
will be prepared at the conclusion of planned work. 
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ABSTRACT 


Transient and steady analysis are used to determine the psychrometric 
conditions that develop in large aboveground fallout shelters ventilated with 
unconditioned ambient air. Thebe analyses consider the shelter size., geometry 
and construction, the psychrometric condition of the ambient weather, and the 
various metabolic and nonmetabolic heat loads to the shelter air'. The results 
of this study indicate that during the hot s.ummer weather, only a small fraction 
of the total energy input.to the shelter is lost through the shelter boundary 3ur 
faces. Thus, the ventilation requirements for large aboveground shelters can 
be obtained by the use of an analysis which neglects the heat loss through the 
shelter boundaries. This means that aboveground shelter ventilation systems 
should be designed to remove the entire thermal load generated within the shelter 

INTRODUCTION 

The fallout shelter must protect the' shelter occupants from radioactive 
fallout and from the detrimental effects of excessive carbon dioxide levels, 
insufficient amount of oxygen, and excessive shelter temperatures and humidities. 
The shelter ventilation system's equipment insures that these additional safe¬ 
guards are maintained. For example, the carbon dioxide level in a shelter can 
be kept at an acceptable value by ventilating the shelter with 3 cfm per occu¬ 
pant. But, as much as 60 efm of outside air could be required to keep the 
shelter temperature and humidity down to tolerable limits during the hot summer 
months. In some situations, ventilating with external ambient air will not 




sufficiently control the shelter effective temperature and air conditioning 
equipment will he required. 

To properly select the ventilation equipment, the shelter designer must 
accurately know the .ventilation rate required. The major parameters in this 
determination are the number of shelter occupants, the shelter's geographical 
location, and the shelter’s size and construction. The number,of occupants 
establishes the level of latent and sensible metabolic heat that is produced 
in the shelter. The geographical location determined that statistical nroba- 
bability and frequency of occurrence of the psychrometric condition ,of the ambient 
weather from which the weather design criteria can be obtained. A knowledge 
of the shelter's physical size and construction features defines the heat 
transfer coefficients which.are required to compute the heat transfer to and " ■ 
from the shelter. When these factors and the required shelter effective temper¬ 
ature, are known and incorporated into a prediction technique for the ventilation 
requirements of the shelter, the designer has a basis upon which he can select 
the ventilation equipment. • ... 

The MED .Division has been engaged in the development of such a prediction 
technique under the Office of Civil Defense Contract ODD-OS- 63 -176, Subtask 
1215A*. The eventual'goal of this program is to formulate a simplified procedure 
for predicting the ventilation requirements of shelters. This paper presents 
the results of several studies which are preliminary steps .in the development 
of this simplified procedure. These studies are based upon 

1. the metabolic energy output of approximately 400 Btu/hr-occupant 
given by the sensible and latent heat expression, see Reference 1, 

•'(Mr. F, V, Alien of O. J'a Directorate of Research was project, monitor. 




2, the shelter habitability criteria represented by the ASHRAE effective 
temperature (an empiracally devised index of the various phychrometric 
conditions that produce similar comfort levels), see Appendix A, and 

3. the derived mathematical shelter model which is applicable to large 

. shelters that are either aboveground or belowground, do or do not have 
boundary surface heat losses, and are in a transient or steady-state 
condition of mass and energy transfer, see Appendix B and Reference 1. 

The mathematical model for the shelter is the most comprehensive analytical 
procedure that could be developed without unduly complicating the computational 
procedure. When used as a transient analysis, the model permits the psychrometric 
State of the shelter air to be computed as a function of time. The analysis 
considers 

1. time-varying inlet air conditions, 

2. time varying energy Inputs to the shelter air from equipment and lights, 

3. time varying solar loads which have been transmitted through windows 
into the shelter (i.e., time varying values of transmitted 3olar 
radiation as obtained from the ASHRAE Guide and Data Book), 

4. metabolic latent and sensible energy loads based upon the instantaneous 
psychrometric state of the shelter, see Reference 1, and 

5. shelter boundary surface heat losses (or gains) based upon a one-dimensional 
heat transfer analysis neglecting corner effects. 

The shelter model is based upon the assumption that 

1. the air within the shelter is so completely mixed that all of the 
shelter air .is at one psychrometric condition, 

2. the convective heat transfer coefficient for each external boundary 
surface of the shelter is not a function of temperature and has a 
constant value over the surface, 


3 








3. the walls and floors which are internal to the shelter volume are at 
the dry-bulb temperature of the shelter air, 

4. the radiative energy interchanges within the shelter can be neglected, 

5* the solar direct and indirect radiative input to the shelter through 

windows and the equipment and lighting heat loads in the shelter can 
be grouped together as one time dependent load factor which is termed 
the thermal load, 

6. the ventilating air exhausted from the shelter is at the psychrometric 
condition of the shelter air, 

7. the effects of condensation on the walls, floors, and ceilings of the 
shelter can be neglected, ' 

8. the thermal-physical properties of the structural materials are not 
temperature dependent, and 

9 . the solar radiation absorption on opaque shelter bondary surfaces 
can be neglected. 

Thus the shelter is idealized as an enclosed volume in which sensible 
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such a model are the conservation of energy and the conservation of the masses 
of dry air and water vapor. 

The primary aim of this paper is to determine the most simple analytical ■ 
model which can predict the psychrometric conditions which develop in above- 
ground shelters. To evaluate this model, actual shelter test data are compared 
to computed data; and variation of parameter approaches are applied to several 
types of shelters to generalize the results. The analytical shelter models 
studied in the paper are 

1. aboveground adiabatic boundary shelter model (no heat transfer 
through the shelter boundary surfaces) with time-dependent parameters 
(l.e., the state of the inlet air and the thermal loads added to the 
shelter air), 

2. aboveground nonadiabatic boundary shelter model with time-dependent 
parameters, 

3. aboveground adiabatic boundary shelter model with steady-state' 
(constant with time) parameters, 

4. aboveground nonadiabatic boundary shelter model with steady-state 

parameters. ■ '? 

RESULTS OF SHELTER ANALYSIS 
Transient Analyses of Shelters 

The time-history of the temperature and humidity within a shelter is pre¬ 
dicted by the transient analysis which considers the time varying mass and 
energy balance about the shelter. The validity of this analysis is established 
by comparison of several sets of analytical computations with data from 
actual shelter tests. The shelter tests chosen for the comparison are the 
MRD Wilmington, North Carolina test #7 (Reference 2), the University of Florida 



Central Stores Building test Phat. . IV (Reference 3 ), and the MRD Houston, Texas 
test II (Reference 4). The Wilmington shelter is a 210-man aboveground shelter, 
the Florida shelter is a 250-man basement shelter, and the Houston shelter is 
a 290-man basement shelter. Using the observed inlet air data and the shelter 
dimensions and construction demils, the transient analysis computed the dry- 
bulb and effective temperatures for the shelter air as. a function fj.f time, see 
Figures 2, 3 , and . 4 . The instantaneous calculated values are generally within 
2°F of the experimental data for the shelter dry-bulb temperature and within 
1°F for the shelter effective temperature. 

When the shelter is assumed to have adiabatic boundaries (no heat transfer 
through the shelter boundary surfaces), the analytical and experimental results 
agree on the 24 hr average to within 2°F for the dry-bulb temperature and within 
1°F for the effective temperature. The adiabatic boundary results are consis¬ 
tently at or greater than the experimental values, This means that all of the 
tested shelters lost energy through their boundary surfaces and that the 
transient analysis siightly overestimates the amount of this energy loss. 
Steady-State. Analysis of Shelters with Bo u ndary Heat Loss 

The time-average psychrometric condition of the shelter can be determined 
by the steady-state analysis which considers the time-average values of the 
psychrometric condition of the inlet air and the heat loads generated within 
the shelter. The heat loss from the boundary surfaces of an aboveground shelter 
is determined by the temperature difference that exists between the shelter air 
and the ambient air external to the shelter and by the heat loss coefficient, UA. 
The UA value for a shelter is determined by the size, geometry, and composition 
of its walls, floor, and ceiling; and its value is determined by 
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Time-Average Temperatures of Shelter 
Environment* 

Dry-Bulb,T db Effective, 

- : Experimental Data 88.8°F 8U.5°F 

- : Adiabatic Boundary 90.5°P 8b.6°F 

Computation 

Nonadiabatic Boundary 88.3°F 84,5°F 

Computation 

*For 20.0 CFM/OCCUPANT Ventilation Rate, 83.9°F Average Inlet Air Dry-Bulb 
Temperature and 80'’F Average Inlet Air Effective Temperature. 

FIGURE 2 - COMPARISON OF COMPUTATIONS AND EXPERIMENTAL 
DATA FOR WILMINGTON, N,C„ TEST 7 
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Time-Average Temperatures of Shelter 
Environment* 

Dry-Bulb, T^ B Effective, 

- : Experimental Data 88.6°F 83.2"F 

- : Adiabatic Boundary 91.5°P 84.2°F 

Computations 

; Nonadiabatic Boundary 87.5 n F 82.0^ 

Computations 

*For 13.9 CFM/OCCUPANT Ventilation Rate, 8l.5"F Average Soil Temperature, 
85°F Average inlet Air Dry-Bulb Temperature and y8.3°F Average Inlet Air 
Effective Temperature 

FIGURE 3 - COMPARISON OF COMPUTATIONS AND EXPERIMENTAL 
DATA FOR CENTRAL STORES BUILDING GAINSVILLE, 

FLORIDA 
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Time-Average Temperatures of Shelter 
Environment* 


Dry-Bulb, T 

--: Experimental Data 88.7°F 

- ! Adiabatic Boundary Computation 90,5°F 

Nonadiabatic Boundary 87.7°F 

Computation 


Effective, T_ 
' E 

84.8°f 

86.1°P 

84.5°F 


*For 12.8 CFM/OCCUPANT Ventilation Rate, 87.9°F Average Soil Temperature, 
83.0°F Average Inlet Air Dry-Bulb Temperature and 79-°°F Average Inlet Air 
Effective Temperature 

FIGURE 4 - COMPARISON OF COMPUTATIONS AND EXPERIMENTAL 
DATA FOR HOUSTON, TEXAS TEST II 
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where: F 

m - number of boundary surfaces 

Uj = overall heat transfer coefficient of boundary surface J* 

A = surface area of boundary surface i 

U 

p = number of shelter occupants 
To evaluate the accuracy of the analysis, the Wilmington aboveground 
shelter test data are compared to the results computed by the steady-state 
analysis with boundary heat loss. The comparison shows that the time-average 
shelter dry-bulb temperature is determined to about 1.5°F and the time-average 
shelter effective temperature is determined to about 2°F, see Table 1, 

TABLE 1 


Ventilation Rate, cfm/occupant 
Average Inlet Dry-bulb Temp,, °F 
Average Inlet Effective Temp., °P 
Average Inlet Relative Humidity 

Average Shelter Dry-bulb Temp,, °F 
Average Shelter Effective Temp., °F 


Values Predicted by Steady-State Analysis with Boundary Heat Loss 


1. ABOVEGROUND SHELTER 

TEST** 



Test Nos. 


5 

7 

8 

10A 

9.0 

20.0 

7.0 

13.0 

83.1 

83.9 

77.5 

90.4 

78.6 

80.7 

74,4 

80.1 

0.68 

0.78 

O .72 

0.4l 

90.0 

■ 88.8 

' 89.8 

94.9 

85.7 

81)..5 

86,1 

85.4 


Shelter Dry-bulb Temp., °F 
Shelter Effective Temp., °F 


*As obtained from ASHRAE Guide and Data Book - Fundamentals and Equipment , 
American Society of Heating, Refrigeration, and Air Conditioning Engineers, Inc,, 
1963 , pps. 391-te8. 

**For 210-man aboveground shelter of Reference 2. 




These results are similar to the agreement obtained with the transient analysis 
If the steady-state shelter values of test #7 are compared to the time-average 
shelter values of the transient study for test #7, the results are within 2°F 
of each other. This shows that the steady-state shelter condition approximates 
the time-average shelter condition determined by the transient analysis. 

The steady-state analysis with boundary heat losses was used to determine 
the psychrometric condition of the shelter as a function of the ventilation . 
rate per occupant, the equipment and lighting load per occupant, the psychro¬ 
metric conditions of the inlet air, and the heat loss coefficient, UA, see 
Figure 5, 6, and 7. The shelter dry-bulb temperature variation due to changes 
in the relative humidity of the inlet ventilation air is negligible (less 
than 0.01°F). This is not true of the effective temperature of the shelter 
and therefore, inlet air relative humidities of 15 and 80 $ are presented. 

The main effect of the heat loss coefficient, UA, is to decrease the 
shelter dry-bulb and effective temperatures as the UA value increases. But, 
the temperature reduction decreases as the flow rate is increased, see Figure 8 
Generally, less than 20$ of the total energy input to the shelter is lost 
through the shelter boundaries; however, this percentage can become as high as 
30-50$ for external ambient dry-bulb temperature below 70°F and flow rates 
below 20 cfm/occupant. 

Steady-State Analysis of Adiabatic Boundary Shelters 

The adiabatic boundary shelter model differs from the nonadiabatic 
boundary shelter model in that it neglects any heat loss through the shelter 
boundary surfaces. This is not an unrealistic assumption, because the heat- 
losses of the aboveground shelter can be a small percentage of the total heat 
input to the shelter. In the belowground shelter, a quasi-adiabatic boundary 
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Neglecting solar absorption on walls. 

Equipment and lighting Load: 

-o mu/at-occuPANT 

- 50 BOT/HR*OCCUPANT 

-100 btu/hr-occupant 


In All These Calculations 
Metabolio Loada*400 BTO/HR-OCCUPANT 
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30 


kO 50 60 


Ventilation Plow Rate, F]_, CFM/OCCUPANT 


FIGURE 5 - NONADIABATIC BOUNDARY SHELTER DRY-BULB TEMPERATURES FOR ALL INLET 
AIR RELATIVE HUMIDITIES AND UA = 10 BTU/HR-“F-OCCUPANT 


Inlet Air Dry-Bulb 













Inlet Air Dry-Bulb 











I 

I 

I 

I 

I 

I 

I 

I 

I 

i 

I 

I 

I 

I 



FIGURE 7 - NONADIABA'TIC BOUNDARY SHELTER EFFECTIVE TEMPERATURES FOR 80f 0 
INLET AIR RELATIVE HUMIDITY AND UA " 20 BTU/HR-°F-OCCUPANT 
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condition is reached when the wall temperature approximates the shelter average 
dry-bulb temperature. For example, consider the Houston basement shelter tests. 
In these shelter tests the time-average soil temperature was 8 t. 9°F whereaB 
the time-average shelter dry-bulb temperatures varied from 82 to 8f°F. When 
the Houston shelter is considered as an adiabatic boundary shelter by .the 
steady-state analysis, the computed data agree with the test data to within 2 6 F 
in dry-bulb temperature and 1,5 °F in effective temperature, sec Table 2. 



TABLE 2 




SUMMARY OF 

HOUSTON SHELTER 

'TEST* 





Test Nos. 


Observed Data 

I 

II 

III 

IV 

Ventilation Rate, cfm/occupant 

9.25 

12.8 

9.25 

18.5 

Average Inlet Dry-bulb.,Temp., °F 

82 

83 

77 

79 

Average Inlet Effective Temp., °F 

78 

79 

76 

77 


Average Inlet Relative Humidity 
Average Shelter Dry-bulb Temp., °F 
Average Shelter Effective Temp., °F 


O .77 0.75 

89.5 88.7 


Average Shelter Effective Temp., °F 87 84.8 

Values Predicted by Steady-State Adiabatic Boundary Analysis 


0.79 

84 

82 


Shelter Dry-bulb Temp., °F 
Shelter Effective Temp., °F 


90.9 90.1 
87.4 85.9 


In test II, the steady-state and transient results agree to within Q.5°F when 
the shelter is considered to have adiabatic boundary walls. This again confirms 
the relationship of the time-average values of the transient analysis with the 
adiabatic boundary steady-state analysis. 

The shelter's dry-bulb and effective temperatures were computed by the 
adiabatic boundary steady-state analysis for the same range of parameters used 
with the nonadiabatic boundary analysis, see Figures 9 > 10, and 11. When 
these data are compared with the data for the nonadiabatic boundary shelter, 
the inlet air conditions that produce high rates of energy loss are found to 


*For 290 -man belowground basement shelter of Reference 4. 














Equipment and Editing Load: 

--0 hrj/hr-occupant 

- 50 btu/hr-occotant 

--“100 bm/hr-occupaht 

In All Tbeae Calculations 
Metabolic I/oad^LOO HIU/HR-OCCUPAMT 



Ventilation Plow Rate, P 1; CFM/OCCUFANT 


ADIABATIC BOUNDARY SHELTER EFFECTIVE TEMPERATURES FOR 
15$ INLET AIR RELATIVE HUMIDITY 









Inlet Air Dry-Bulb Temperature 














occur when energy loss is the least needed to control the shelter environment. 

For example, at inlet dry-bulb temperatures below TO°F, the adiabatic boundary 
shelter's effective temperature never exceeds 85 °F except at flow rates of 6 cfm/ 
occupant and below, see'Figures.10 and 11. 

If the effective temperature index is accepted as the shelter habi- 
ability criteria, a psychrometric chart can be used as a means of determining 
the relationship between the psychrometric condition of the inlet air and the 
ventilation flow rate, Bee Reference 5- For instance, if an average effec¬ 
tive temperature of 8 5 0 F is chosen as the habitability limit, the loci of 
the average psychrometric conditions of the inlet air that produce an 85 6 F 
effective temperature in the shelter can be plotted for various ventilation 
rates, see Figure 12. 

Equipment and Lighting Load = Zero 
Metabolic Load ® 400 Btu/hr-occupant 
Steady-State Adiabatic Boundary 
Shelter Model 



FIGURE 1.2 - VENTILATION RATES REQUIRED TO MAINTAIN AN AVERAGE EFFECTIVE 

TEMPERATURE OF 85°F 




Local Ventilation Rates 


In order to determine the ventilation rate required for any giyen shelter, 
weather design criteria must be selected. Several weather data studies are 
currently available upon which the weather design criteria could possibly be 
based (e.g., see References 6, 7, and 8), but the weather design criteria to 
be selected for a shelter must be based upon a detailed study of the relation¬ 
ship between the psychrometric state of the ambient weather and the effect that 
ventilation air from this ambient has on the shelter environment when actual 
weather data are considered. This relationship will be influenced by the Interval 
of time chosen for the comparison. Possibly none of the available weather data 
studies will be applicable and a new set of criteria will have to be established. 
l‘he MRD Division is presently engaged in a program for the Office of Civil 
Defense to study the relationship of periods of actual weather on shelter envi¬ 
ronment. In addition, more information must be gathered concerning the effect 
of environmental conditions on the human body; particularly, when these con¬ 
ditions are changing with time and the body is under a high level of emotional 
stress. However, this paper has shown that the regional ventilation requirements 
for shelters can be determined once the weather design criteria and habitability 
criteria are known, 

CONCLUSIONS 

The transient analysis is able to predict an aboveground shelter's instan¬ 
taneous dry-bulb temperature to within 2°F and the shelter's instantaneous 
effective temperature to within l"F, The time-average values of a shelter's 
dry-bulb and effective temperature are predicted by a steady-state analysis to 
within 2°F. Either of these analyses can be used to determine the environ¬ 
mental condition of a shelter with a reasonably high degree of accuracy. The 
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steady-state analysis has established that the energy loss that can occur 
through the boundary surfaces of an aboveground shelter is generally less than 
20 $ of the total thermal energy introduced into the shelter air during hot 
summer weather. All of these results indicate that the.mechanism of heat loss 
through the boundary surfaces of an aboveground shelter cannot be depended 
upon to remove energy from the shelter during hot weather. At most, boundary 
surface heat loss should be regarded as a possible safety-factor in ventilation 
system design. It is therefore recommended that the ventilation systems for 
aboveground shelters be designed to remove the entire thermal load generated 
within the shelter. , 

The flow rate predicted by the analysis of a shelter without boundary 
surface heat loss is the maximum flow rate that can be required (assuming the 
solar radiation absroption effects on opaque shelter boundary surfaces are 
negligible *)5 and thus provides a means for establishing an upper limit on the 
sise of the ventilation equipment for an aboveground shelter. The reasonable 
agreement between the calculations for a shelter without boundary surface heat 
loss and the shelter test data insures that the equipment size based upon this 
ventilation rate is not overly conservative. 

The paper has shown that with metabolic head load data and weather design 
criteria the ventilation requirements for aboveground shelters (e.g,, the 
shelters surveyed in the National Fallout Shelter Survey) can be determined 
analytically. The reliability of these predictions is primarily dependent upon 
the reliability of the metabolic and weather design data used in the calculation 
procedure. 

*This assumption is presently under detailed evaluation by the MRD Division. 
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DERIVATION OF EFFECTIVE TEMPERATURE EQUATION 


Many sets of psyehrometric conditions provide the human body with similar 
levels of physiological comfort. The effective temperature scale has been 
established as an empirically derived index of the various psyehrometric condi¬ 
tions which produce similar eomfort levels. The effective temperature is a 
function of the dry-bulb temperature, wet-bulb temperature, and velocity of the 
air in which a person resides. Generally, these data have been presented in 
nomographic form or by approximating equations. For the purposes of this study, 
the nomograms were inconvenient to use since they were not in equation form, 
and the existent equations were not of sufficient accuracy. As a result, an 
accurate mathematical expression for effective temperature was derived. 

The derivation was based upon the ASHRAE nomogram* for effective temperature, 
see Figure Al. Because the ventilating air velocities in shelters is generally 
low, only the 20 ft/min flow velocity curve of the nomogram was considered. 

This eliminated the a.ir velocity as a parameter in the determination of effec¬ 
tive temperature. A schematic of the curve AA which represents the 20 ft/min 
flow velocity in the nomogram is shown in Figure A2. Curve AA is linearized 
by a straight .line which intersects the wet-bulb (T, m ) and dry-bulb ('I') temper- 
ature scales at the points M and N respectively. By geometric similarity 


P-Q M-Q 

N-0 ' M-0 


(Al) 


and 


P-Q R-Q 

S-T R-T 


(A2) 


* AGHRAE Guide and Data Book - Fundamentals and Equipment , American Society of 
Heating, Refrigeration, and Air Conditioning Engineers, Inc., 190.3, p. Ill 


A2 















Substituting the above values into Equation (A3) yields 


EFF 


107 (I DB 

62.3 + (T 


• V * 62 - 3 <W 


DB 




(A4) 


This relationship is limited to low air velocities and is restricted to the 


temperature range of 

'.5-P < T M < 110-F 


and 

30°F < T wb < 100°F 

The temperature range restrictions are imposed due to the fact that the T-,™ 

ii). 1 i 

curve of the nomogram cannot be considered a straight line beyond these temper¬ 
ature limits. Table A1 shows that Equation (A4) has an average error of less 


than one percent. 

TABLE A1 ERROR IN CALCULATED ERRECTIVE TEMPERATURES BY EQUATION (A4) 


T . °F 
DB’ 

V ° F 

T eff’ Cai ' d by 
Equation (A4) 

T eff’ read from 
Figure (Al) 

Percent Error 

110.0 

98.0 

99.4 

99-0 

0.45 

110.0 

90.0 

94.2 

94.0 

0.21 

110.0 

30.0 

73.5 

74.0 

' 0.68 

90.0 

80.0 

03.9 

83.9 

0 

90.0 

60.0 

75-5 

75.6 

0.13 

90.0 

4 o.o 

70.0 

70.0 

0 

80.0 

60.0 

71.5 

71.5 

0 

80.0 

4 o.o 

66.3 

66.3 

0 

80.0 

30.0 

64.7 

64.3 

0.62 

70.0 

60.0 

66.6 

66 .5 

0.15 

70.0 

50.0 

63.9 

64.0 

0.15 

1 + 5.0 

30.0 

45.1 

45.0 

0.22 


The effective temperature values are correlated to the relative strain 
indicies with a 0.5 relative strain index being essentially an 85 "F 
effective temperature, see Figure A3- Almost all people are comfortable at 
zero relative strain index and physical failure is rapid and severe at 1.0 


A 3 





imTii 












relative strain index.* This comparison illustrates the convertability of the 
effective temperature index to other indices of physiological comfort. 


*11. 11. Lee and A. llensehel "Evaluation of Thermal Environment in Shelters", U.S. 
Department of Health, Education, and Welfare, TR-8, August, 1963 
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APPENDIX B 


ANALYTICAL MODEL 



B.l DERIVATION OF ANALYTICAL SHELTER MODEL 


The analytical model for a shelter is formulated with the philosophy that 
whether a shelter is above or below ground and whether it is single story or 
multi-story, the shelter can be represented by a single mathematical model. 
Furthermore, the model is developed with the idea that it should be the simplest 
model that can be devised and still adequately explain the phenomena. In accord 
with this philosophy, the assumptions are made that 

1 . the air within the shelter is so completely mixed that all 
of the shelter air is at one psychrometric condition, 

2 . the convective heat transfer coefficient for each external 
boundary surface of the shelter is not a function of tempera¬ 
ture and has a constant value over the surface, 

3 . the walls and floors which are internal to the shelter volume 
are at the dry-bulb temperature of the shelter air, 

4. the radiative energy interchanges within the shelter can be 
neglected, 

5 . the solar direct and. indirect radiative input to the shelter 
through windows and the equipment and lighting heat loads in the 
shelter can be grouped together as one time dependent load factor 
which is termed the thermal load, 

6 . the ventilating air exhausted from the shelter is at the 
psychrometric condition of the shelter air, 

7. the effects of condensation on the walls, floors, and ceilings 
of the shelter can be neglected, and 

8 . the thermal-physical properties of the structural materials 
are nob temperature dependent. 


it’ 



The model of the shelter is established by these assumptions. Thus the shelter 
is idealized as an enclosed volume in which sensible heat, latent heat, solar 
heat, equipment heat, and ventilating air are introduced'and from which air is 
exhausted and energy is lost. The governing principles for such a model are 
the conservation of energy and. the conservation of the masses of dry air and 
water vapor. 


The conservation of energy requires that 


the rate of change of enthalpy the rate of c 

in the ventilating air + in the shelter air 


J 


the rate of metabolic heat 
input by the shelter occupants 


j ["the rat 

iS + I load to 

and integrating Equation (Bl) over the time increment At 


the rate of change of enthalpy - ] 


ate of thermal”] 
to the shelter 


' 


the rate of energy transfer across" 
the shelter boundary surfaces 


(Bl) 


(“l - V • V> * S. ♦ «T - % - 0 < K > 

whore: 

- enthalpy of ventilating air entering the shelter, Btu 

Hg = enthalpy of ventilating air leaving the shelter, Btu 

H„ . = enthalpy of shelter air at beginning of time increim Ai , if o 

b, 1 

H 0 „ = enthalpy of shelter air at end of time increment A".' , btu 
b, d 

- energy due to metabolism of shelter occiipants, Btu 

~ energy input due to thermal load (see assumption #5.> p. B2), Btu 

- energy loss or gain through shelter boundary surfaces, Blu 
a 



The conservation of the mass of water vapor necessitates that 


t 

\ 


the rate of change of water - ] 


vapor in the ventilating nirj + j^vapor 1 


Ptbe rate of chang 


in the shelt 


;e of water 
Lter air 


the rate of water vapor 

introduced by the shelter oceupantsj = 0 


s 


and integrating Equation (B3) over 'the time increment At 


( Vl ' M V,2 ) + ( \s,l ' ^,S,2 ) + *V/] 


(B 3 ) 


(bU) 


where: 


My j -• mass of water vapor in the entering ventilating air, lb 

My g = mass of water vapor in the exhausting ventilating air, lb 

My g ^ = mass of water vapor in the shelter at the beginning of 

’ ’ the time increment kr , lb 

My p - mass of water vapor in the shelter at the end of time 

’ ’ increment A T , lb 

My M « mass of water vapor introduced by the shelter occupants, lb 
The conservation of the mass of dry air demands that 


the ra 
air in 


rate of change of dry 
the ventilating air 


y the rate of change of dry 

ir + air in the shelter air •- 


and integrating over the time increment Ar 


(M . - M ) + (M - M _ c 

a,l a,2 a,S,l a,S,2 


) = 0 


0 

(B 5 ) 

(b6 ) 


where: 


a, 1 
^a, 2 


M 


a, S, 1 


a, o, 2 


mass of dry air in the entering shelter air, lb 

mass of dry air in the exhausting shelter air, lb 

mass of dry air in the shelter at the beginning of the 
time increment, lb 

mass of dry air in the shelter at the end of the time 
increment, lb 




Denoting 


AQ = li ± + H Sjl + Qn + - Q b (B7) 

Vo = Vl + V S, 1 + V M (b8 ) 

M a,0 = M a,l + M a,S,l (B9) 

and substituting Equation (B7) through (B9) into Equations (B2), (b 4) and 
(B 6 ) gives. 

, AQ = Hg^g + Hg (BIO) 

^,0 " \2 + \S;2 (Bll) 

M a,0 = M a,2 + M a,S,2 (B12) 

wherei 

H S,2 = M a,S,2 h a,S,2 + Vs,2\s,2 (B13) 

, H 2 = M a,2 h a,2 + V2 h V,2 (Bl4) 

From the assumption that the exhaust air has the same psychrometric condition 
as the shelter air 

h a,2 = h a,S,2 (B15) 

and 

h V,2 \s,2 (Bl 6 ) 

With Equations (Bll) to (B1.6) substituted into Equation (BIO) 

= M a,0 h a,S,2 + Vo h V,S,2 (B17) 

The specific enthalpies of dry air and water vapor are given by 

h = 0.24 (T nn + 459-69) (B18) 

3, 19 Jj 

and 

h - 1061 + 0.444 (T„J (B19) 

v .DB 

for 32°F < T db < 150°F 

b:> 



where: 


= dry-bulb temperature, °F 

Combining Equations (Bl?) to (B19) and solving for the dry-bulb temperature 
results in 


Vo 

DB,2 0.444 My 0 + 0.24 M & ~ 


(B20) 


where: 


dry-bulb temperature of shelter air at end of time increment 

At , y 

The following relationships are evident: 




. M a,2 = 

6 ° P 2 p a,2 AT 

(B21) 


^V, 2 

60 *2 p V, 2 AT 

(B22) 


M a,S,2 

V P a,S,2 

(B23) 

and 

\s,2 = 

V P V,S,2 

(B24) 


where: 


V = shelter volume, ft 
At = time increment, hr 

O 

0 = density of dry air leaving shelter, lb/ft J 

Qi) cL 

,Oy g= density of water vapor leaving shelter, lb/ft 3 


p 0 - density of dry air in shelter at end of time increment, 

a ’°’ 2 .lb/ft 3 


P 


V,S,2 


t 


2 


density of water vapor in shelter at end of time increment 
lb/ft 3 


volumetric flow rate of exhaust air, 


ft 3 /min 


but by assumption 



/ a, 2 

p a, S, 2 

(B25) 

and 

P V,2 

' P V . 0 , r> 

(B26) 



Bo 







Substituting Equation (B2l) through (B26) into Equation (Bll) and (B12) yields 

\,0 " P V,2 (6 ° F 2 Ar + V) (B2Y) 


and 


".,0 ■ <“ V * V > 

When Equations (B27) and (B28) are solved for Eg, they result in 

M - p V 
a, 0 'a, 2 


5oSr7 a>2 

ny .o ~ p v,2 v 


6oat p. 


(B28) 

(B29) 

(B30) 


V,2 


The densities of the dry air and water vapor are assumed to obey the perfect 
gas law. Thus; 

P 


V 


and 


where: 


v ~ R v (T db + 459 . 89 ) 


p - p 
B V 

P a R a (T DB + 1+59.69) 


(■B31), 


(B32) 


P B = barometric pressure, lb/ft 

2 

Py = partial pressure of water vapor, lb/ft 

Ry = gas constant for water vapor, ft-lb/lb mole -°R 

R = gas constant for dry air. ft-lb/lb , -°R 

a ' mole 

The partial pressure of water vapor is expressed by 


r P„ 


(B33) 


where: 


relative humidity 


P_ ■= saturation pressure of water vapor, lb/l't e 


P g -- '3. .132e° ’( t bb ) 


(B3 j i: 


fur 32“f< T np 


1'30°E 


BY 




which is an analytical expression for the tabulated values of saturated vapor 
pressures as a function of dry-bulb temperature (Ref. 8 ). 

From Equation (B30) through (B 33 ) with 


and 


»« * 53 - 35 n - lb / 1 W R 


”v - 85.71 


the following expressions can be obtained: 

P V 


P„ a r [5.132e°‘ 329 ^ T DB^] 




and 


P, 


l 


a (,t db + d>y.oy; • B 

Equating Equation (B27) and. (B28) gives 


M. 


a, 0 


“v.O 


(B35) 


[ 0.05987 r e°' 03 i! 9 < V] (B36) 

[O.Ol8?44P - 0.0962 r e °•°329(T DB )] ( 537 ) 


p a ,2 PV ,2 

Evaluating Equations (B36) and (B37) at the shelter dry-bulb temperature, 
[t db g J,and substituting the results into Equation (B 38 ) gives the partial 
pressure of the water vapor [Py].as 


(B3B) 


(B39) 


v,o 


And from Equation (B33)> the relative humidity of the shelter at the end of 


the time increment is given by 


r - 


(B4o) 


Substituting Equations (B34) and (B39) into Equation (b4o) results in 


Bo 





V 


(B4l) 


5-132e (0.0329(T db )> [l * 0 . 6221 , i°l 

1 “v,o J 


With these equations, the psychrometric conditions within the shelter can be 
determined as a function of time through the use of the following procedure. 

1. Determine the psychrometric conditions of the air 
introduced into the shelter and the air within the • 

, shelter along with the heat inputs and losses of 

the shelter for the time interval At 

2. Compute the quantities AQ, M,, Q , and M a Q from Equations 
(B7)> (b 8), and (E9) respectiv4ly. ’ 

3. Assume that these quantities do not change over the time 
interval AT, and compute the dry-bulb temperature of the 
shelter.at the end of the time interval from Equation (B20). 

4. Calculate the relative humidity of the shelter air at the 

end of the time interval from Equation (B4l). \ . 

5. Set the psychrometric conditions of the shelter at the end of 

the time increment equal to those in the shelter for the beginning 
of the next increment, and repeat the entire procedure. 

By the continuous application of this procedure, the shelter's psychrometric 

condition can be obtained as a function of time for any time period. The 

computation method that has been developed constitutes a transient analysis 

of the shelter environment under the influence of time varying parameters. 


The shelter's psychrometric condition is described in the analysis by its 
dry-bulb temperature and relative humidity. Instead, the dry-bulb and wet-bulb 


or effective temperatures can be used. The wet-bubl temperature is determined 


from the Carrier equation. 


l ' P B " P s'^ T DB ~ 

2800 . 0 - 1.3 (t wb ) 


(b42) 


J. H. Carpenter, "Fundamentals of Psychrometries", 
Technical Note T 200-20, Syracuse, New York, 1962 , 


Carrier Corporation, 
p. 12. 
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where: 


2 

Pg 1 = saturation pressure at wet-bulb temperature, lb/ft 

T,_ = wet-bulb temperature, °F 

Wi3 

with P 1 evaluated by Equation (B34) 

P g ' = 5-132 e 0,0329 ^WB) (BU-3) 


and obtained from Equation (B39)• The results is a transcendental equation 
that can be solved for the wet-bulb temperature, Tyg. When the value of the 
wet-bulb temperature is known, the effective temperature can be determined 


from equation (A4). 


10 7.5 < t db - W + & ' 3 <V 


= 2 ' 3 * < t db - W 


As mentioned in Appendix A, this relationship is limited to low air velocities 
and is restricted to the temperature range of 

45 e F < T' < 110 n F 

— JJij 

and 

30 °F <T ra < 100 °F 


The heat loss or gain of the shelter boundary, Q^, over a time interval At 
is determined by the temperature of the inner surface of the boundary, T^. In 
the computations, the value of T^ is obtained from the temperature distribution 
that existed through the boundary during the previous time increment. This 
temperature distribution is deduced from a transient analysis of the energy 
transfer in the boundary. In order t.o accomplish this analysis, the boundary 
is divided into a finite number of slabs of thickness AX, except for the inner 
and outer surfaces which are made into slabs of thickness A X/2, with each slab 
assumed to be at a single temperature. This temperature is assigned to the 
midpoint of each slab, except for the innermost and outermost slabs which are 








assigned temperatures at the external surfaces of the slab (see Figure B.l). 
These locations of slab temperature are termed nodal points. The external 
nodal points of the boundary transfer energy by'convection with surrounding air 
arid by conduction with the next internal nodal poipt. All of the rest of. the 
nodal points transfer energy by conduction with the two adjacent nodal points. 
The energy balance about an- internal nodal point in is 


"energy conducted"" 


"energy conducted to 


[""energy stored in ] 

to point m from 

... - 

point (m+1) to 

-•.. 

[the slab of point ml 

point (m-l) 


point m 


(34*3) 


(T -T J 
v m m-1' 


(T ^-T ) 
m+1 nr 


AX>C (T '-T ) 
pm m 


and .the temperature at the nodal point is 



(B47) 

(b48) 



k = thermal conductivity of composite boundary material, 
Btu/hr-ft-°F 

y ' o 

p - density of composite boundary material, lb/ft J 

C - specific heat of composite boundary material, Btu/lb-°F 

Ax - thickness of boundary slab, It 

T - temperature of nodal point in at beginning of time 
111 increment, "F 

T ^ = temperature of nodal point m-.'l at beginning of time 
increment, ‘‘F 










T = temperature of nodal point m+1 at beginning of time 
m increment, °F 

T 1 = temperature of nodal point m at end of time increment, °F 

The energy balance about the inner surface of the boundary is 


energy convected 
from shelter air 
to boundary 


energy conducted to 
nodal point 2 from 
_nodal point 1 


energy stored in 

the slab at 

nodal point 1 (b49) 


' (t 2 -t ) 

h i A ( T DB,2 * T l> + 


f- A ^ (T l' ' T l ) 


and the temperature at the inner boundary surface is 


T i = ■ Wl»,2 + B 3 T 2 + (1 'V B 2 B 3 )T 1 


with B^Bg+l) < 1 


where: 


h i A x 

k 


B„ = 2B, 


h = heat transfer coefficient on the inner surface of the . 

1 boundary, Btu/hr-ft 2 -°F. 

T,,,, -= temperature of shelter at beginning of time increment, °F 

JUBj d 

= temperature of nodal point at inner surface of boundary at 
beginning of time increment, °F 

T = temperature of nodal point 2 at beginning of time increment, 

2 o F 

' = temperature of nodal point at inner surface of boundary at 
end of time increment, °F 

The energy balance about the outer surface of the boundary n is 


energy conducted 
from (n-l) nodal 
point to outer 
surface point n 


na.1 ambient air 


energy convected energy absorbed energy stored 

from outer surface from solar radi- in the slab 
point n to exter- + ation in the slab -- at the outer 
nal ambient air at the outer ^surface 

surface 


surface 


BJ.3 



or 


- T -~— - - ‘.‘'WWb 


AX 


and the temperature at the outer boundary surface is 

T ' 


where : 


At 

B, 


Vn-l * V 3 \ * ^'h'hhK * % 


With B 3 (B 4 + l) < 1 


B, = 


h 0 A * 


(B54) 

(B55) 

(B56) 


h 

o 


n-1 


= heat transfer coefficient on the outer.surface of the 
boundary, Btu/hr-ft- -°F 

=• intensity of solar radiation for shelter location and 
given calender date per unit area, Btu/hr-ft^ 

= temperature of nodal point n at outer surface of boundary 
at beginning of time increment, °F 

temperature of nodal point n-1 at beginning of time 
increment, °F 


T. " = temperature of.external ambient air at beginning of time 
increment, °F 

T ' = temperature of nodal point n at outer surface of boundary ■ 
at end of time increment, °F 

f = absorptance of outer slab surface for solar radiation 

The boundary surface heat loss as defined by Equations (b 47), (B5l), and 
(B55) is applicable to the one dimensional heat transfer from a boundary sur¬ 
face. If the boundary surfaces of a shelter are large, corner effects are 
negligible and each boundary surface ce.n be considered to be conducting energy 


from or to the shelter uni-directionally; However, the problem in solving for 
the heat energy loss or gain from the shelter boundary surfaces is dependent 
upon the designation of a temperature in Equation (B55) and the properties 
of the boundary surface, i.e., density (p ), thermal conductivity (k), specific 



heat (C ), and thickness ( X) with Av = —— 
v p" ' ' A n-1 


where n = number of nodal points. 


To simplify the analysis, the shelter boundary surfaces are grouped together 
to form one slab surface exposed to the temperature of the shelter on one side, 
and the temperature on the other side. The procedure of replacing all of 
the separate boundary surfaces by one slab surface is based upon the determination 
of area weighted average values for the various properties of the boundary sur¬ 
face. The property values of each boundary surface are based upon thickness 
weighted average property values of the materials that compose each of the 
boundary surfaces. For example, for any boundary surface property X, the area 

weighted average value of' property X for the slab surface is given by 

r(j) m A 

X - [ 2 2 { ^L) (X pXp ) ] / A 


P-1 




(B57) 


where 2 


A. 

J 


- surface area of boundary surface ,j 
m = number of boundary surfaces 

Xp = property X of the material p in the boundary surface j 
X p = thickness of the material p in the boundary surface j 
r(,-)) = number of materials in boundary surface j 

r(j) 


X 


J 


A = 


2 

P“1 

m 

2 

j-1 


X 


so that 


r (j) 

[ S 

P-1 

r(j) 

[ s 


A. 


2 

J=1 

m 

2 

j=l 


A O 1/ i 


P P 


( X p k p )] /A 


(B58) 


(B59) 


B15 



m 


r(j) m 

[ 2 2 

p=l j=l 


(- jI ) ( A ( C ) ,) ] / A 
X • p*• p J J J ' 

J 


(BuO) 


and 


[2 A \ , ]/ A 

J --1 J J ■ 


(B6l) 


The values p , k, C '> and X determine the property values of the composite 
boundary surface slab. 


In general^ the exterior sides of shelter boundary surfaces are exposed 

to three types of environments. The boundary surface can be exposed to the 

ambient weather (e.g., an outside wall), to the soil (e.g., an underground 

wall, the floor), or to a space interior to the structure but exterior to the 

shelter (e.g., a first floor shelter ceiling exposed to a second floor space 

of the building). In each of these situations, the value of T. for fitch 

A 

boundary surface will be different. In general, the area weighted average 

value of T, would be . 

A » 


where: 



2 a;(t.). 

d A'j 

T ~ --- 

A m 

2 A. 


temperature to which the exterior side of boundary 
surface j is exposed 


(B62) 


Letting 


P, 




m 

2 A. 
3 


with j - 1,2,..., 


jl 


(B63) 


Blh 


m 



then 


T 


A 


m 

2 

J«1 




(B64) 


With the fact that the exterior side of the boundary surface can be exposed 
to only three types of environments, then 

*A = Wo + ? 3 ( Vs + Wi vith 1 = ,(B6 5 ) 

where:' 

P q = percentage of shelter boundary surfaces exposed to 
the ambient weather temperature (T^) q 

p = percentage of shelter boundary surfaces exposed to 
solid at the temperature (T^)g 

= percentage of shelter boundary surfaces exposed to 
each of the interior space temperatures (T^)^ for 
each of the q interior spaces. 


However, experience and experimental test results have shown that the 


temperature values of the soil and interior areas adjacent to the shelter can 

-11 •» . 

be closely approximated by 


and 


(Vs 

(Vi 


T s + t se 


(Vo + t se 


(b66) 

(b67) 


where: 


T c - well.water temperature at location of interest as given by Collins.* 

Tgg - shelter dry-bulb temperature at time for which boundary surface 
heat loss or gain is computed. 

Substituting Equations (b 66 ) and (B 67 ) into Equation (B 65 ) gives 


^A ' = P o'Vo f 2 [ T S + Vi] h 2 ['Vo + T Se] 


(B68) 


-K 


W. D. Collins, "Temperature of Water Available for Industrial Use in the 
United States", U.S. Geological Survey Water Supply Paper 520-F, .1925. 





ty-r 

®e 


I 


i 


or 


..vhere: 


T = f-p + — ) (t ) + 

A 2' ''ho 


(p s + Pj) 


T SE + 2 P S T S 


(B69) 


Pj = S P^> total percentage of shelter boundary surfaces 
, 1=1 . exposed to interior spaces 

The value of f^ of Equation (B 69 ) is used as the value of in Equation (B55). 
With these area weighted average property values and the relationship for 
the value of the boundary surface heat loss or gain, Qg, can be computed. 

The analysis also considers the solar radiation that enters the shelter 
through windows and the solar radiation that is absorbed on the exterior surfaces 
of the shelter r The solar energy that comes through the windows is treated as a 
time varying thermal load in the term of Eq. (B2). The magnitude of this load 
is determined by ■' | 

<V«inta,~ * ■%>! i} W°> 


where ■, 


(Q g ) w i n dow = ener Sy transmitted into shelter at time of consideration 
« solar intensity on window j at time of consideration 






7 


J 


= surface area of window j 

= transmittance of window j for solar radiation 
The solar energy that is absorbed by the walls is neglected as an assumption 
in the study of the shelters mentioned in the main body of this report. Therefore 
in Eq. (B5*0, i/j is taken to be equal to zero. Generally, this may be handled in 
the following way. 

The solar energy that is absorbed by the wall is defined by the solar intensity 
Q g times the absorptance of the exterior surfaces for solar radiation, ^ [ see Eq. 
(B54) ]. The value of that is used in Eq. (B5^-) is an area-weighted average 

value of the solar intensities that occur upon each of the boundary surfaces of 


BJ.O 






the shelter exposed to sunlight. Consequently, this solar energy term is a 
function of time and is defined as 

i w/j 


Z A. 


where: 


A. = area of the bondary surface J, ft 
J 

n = number of bondary surfaces exposed to solar radiation 
at time for which § g is required 

(Q ) = intensity of solar radiation incident upon the bondary 2 

s ** surface, J, at time for which Q is required, Btu/hr-ft 


The value of (Q ) is defined by 
8 J 


«.), - «, 


direct ^ j + ^diffuse ^J + ^reflective 


where: 


(Q,, . ) . = energy received by surface j directly from sun, 

direct J Btu/hr/ft 2 

^diffuse = energy from atmospheric scattering which 

J reaches surface J, from all directions in 

sky, Btu/hr-ft^ 

(^reflective K = ener sy reflected from objects facing the 
J surface jy that reach the surface J, 

Btu/hr- fir' 

B. 2 COMPUTATIONAL CRITERIA 


The time increment used in a computation cannot be arbitrarily chosen with¬ 
out introducing instability into the computations. Practice has shown that if 
the volume of the ventilating air introduced into the shelter during a time 
increment is less than a fixed percentage of the internal shelter volume, the 
calculations will be stable. That is, 
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sfeg Ti nflfcji r.sESi ■rx.-TsxxiKSXt* *s 


60 (F ) (At) < C V 


where: 


fe(F 1 ) 


C = a constant coefficient (= 0.1, by experience) . 

At = time increment, hr 

F^ = volumetric flow rate of entering ventilation air, ftV min 

■5 

V = shelter volume, ft 


When a shelter's stability conditions are computed, the above stability 
criterion is the relationship that requires the smallest time increment, At, 
and not the relationships of Equations (B48), (B|52), or (B5^) * Therefore, the 
Bize of the time increment is determined by Equation (B7^)* 

\ ■ ^ <' 

Two other computational procedures that are followed a\/the mathematical 
model of the shelter are 

1) the input data are linearly interpolated for time increments that 

,, are smaller than those for which the data are given, and 

2 ) the relative humidity; of Ijhe shelter is always kept less than or '■ 
equal to unity. If the computed relative humidity is greater than 
unity, the latent heat input to the shelter is reduced until the 
relative humidity is equal to unity, 

B3 - STEADY-STATE SHELTER ANALYSIS 

In the steady-state analysis of the shelter environment, none of the para¬ 
meters vary with time including the psychrometric condition of the shelter. 

This reduces the complexity of the computation program. The transient calcu¬ 
lation procedure can be used to obtain steady-state environmental results if all 


B20 







of the parameters are made constant with tim e ' T! process is an iteration 
procedure which converges to the steady-8tat« Values. The steady-s-ta^ ca i cu . 
latlon Is otherwise identical to Use tran#i@nt shelter model. 

Under steady-state conditions the conservation of energy equation becomes 

(H 1 ' V * % + VS* C ‘ (B75) 

the conservation of the mass of water vapor is 

(M V,1 * Vs** Vm.* 0 • (876) 

and the conservation of dry air mass Le . 


Introducing the quantities AQ,* M^q o 

AQ - H x + Q m + Q t - Q 0 

Vo = Vl + Vm 



H g AQ 



(B77) 

(B78) 

(B79) 

(B 80 ) 

(B8l) 


V,2 * Vo v. , . • . > (B82) 

•\,2 = Vo - - (B83) 

All of the other relationships remain the same as in the derivation of the 
transient analysis.. The shelter's dry-bulb temperature is still define by 
Equation (B20) and, its relative humidity is still defined by Equation (plf-i) with 
the exception that the terms AQ, My q > and M a ^ 0 must be defined as in Equations 
(B 78 ) through (b 8 o). One consequence of the steady-state derivation is 
shelter deminsions and volume do not enter into the calculations of "the shelter's 
psychrometric condition. 
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